Plasmonic nature of van der Waals forces between nanoparticles 
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We propose a new approach to calculate van der Waals forces between nanoparticles where the 
van der Waals energy can be reduced to the energy of elementary surface plasmon oscillations 
in nanoparticles. The general theory is applied to describe the interaction between 2 metallic 
nanoparticles and between a nanoparticle and a perfectly conducting plane. Our results could be 
used to prove experimentally the existence of plasmonic molecules [H'lH and to elaborate new 
control mechanisms for the adherence of nanoparticles between each other or onto surfaces. 

PACS numbers: 42.50.Pq; 73.20.Mf; 85.85.+J; 78.67.Bf 



In recent years, metallic nanoparticles have regained 
new and vivid interest. Especially noble- metal nanopar- 
ticles are used in biological sensing applications Q , imag- 
ing and as optical materials [4|. The optical response of 
nanoparticles is governed by the excitation of their sur- 
face plasmon resonances. Surface plasmons play an im- 
portant role in many fields of physics, for example in the 
enhancement of the transmission of light through metallic 
structures @-@], or for dispersion forces in 2-dimensional 
Wigner-like crystals 0. The properties of surface plas- 
mons and, thus, the optical response of the nanoparticle 
depends critically on their material, shape, and size Q. 
For example, optical forces between silver nanoparticle 
aggregates are known to be enhanced by exciting conve- 
niently surface plasmon resonances [9[. Here we investi- 
gate a similar phenomenon, that is the tailoring of van 
der Waals interaction between nanoparticles via surface 
plasmons. This interaction is mediated not by a photon 
field but by the electromagnetic vacuum fluctuations, and 
is particularly interesting for the design of control mecha- 
nisms of the adherence of nanoparticles among each other 
or onto a surface. 

To describe the van der Waals interactions between a 
nanoparticle and a surface, usually a dipole approxima- 
tion is used [Hl-dHj which is valid only for large enough 
distances between the particle and the surface. Mitchcl 
and Ninham [l^] suggested to use bispherical coordinates 
for the calculation of van der Waals forces between two 
spherical particles. Recently van der Waals forces be- 
tween a spherical particle and planar substrate have been 
considered with the help of expansion of electromagnetic 
fields over usual spherical harmonics while in [Fu l- 
fill the van der Waals energy between realistic metallic 
surfaces was shown to be dominated by surface plasmon 
oscillations for small distances. 

In the present paper we study the van der Waals energy 



between two nanospheres and between a nanosphere and 
a perfectly conducting plane as the energy of elementary 
surface plasmon oscillations in such a system. Plasmonic 
oscillations in this system have a very complicated spatial 
and spectral structure, which can be described with the 
paradigm of plasmonic atoms and plasmonic molecules 
We will show that unbound plasmonic atoms and 
plasmonic molecules give contributions to the van der 
Waals energy which are approximately equal in ampli- 
tude but opposite in sign. As a result, the van der Waals 
energy turns out to be very sensitive to geometry and 
materials of the interface and the nanoparticles. We will 
consider only short distance limit when all spatial scales 
are smaller than the plasma wavelength X p i = to p i/c, 
where lo p i is the bulk plasma frequency. We will also 
assume a local description of matter and use throughout 
the paper the dissipationless form of Drude's dispersion 

law, e (ui) = 1 — the plasma model. 

To find the plasmon spectrum within the quasistatic 
approximation we solve the Laplace equation for poten- 
tials A<& = with usual boundary conditions of continu- 
ity for the tangential component of E = — V<& and the 
normal component of D = eE and find resonant values of 
dielectric permittivity, e„ , n=l,2,3 .., or the correspond- 
ing frequencies. It is obvious that for a specific matter 
with e = e (uj) the eigen-oscillations will occur for u) n 
such that e (w„) = e„. For the plasma model, the eigen- 
frequencies and energy of plasmon fluctuations and the 
van der Waals energy will be equal to 
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The eigen- values of e depend on geometrical parameters 
which allows to find the van der Waals forces by differen- 
tiation over those. Therefore, the problem is to find uni- 
versal eigen-values of the dielectric constant and, then, 
the frequencies of the plasmonic resonances for specific 
materials. We then apply this approach to the van der 
Waals interaction between 2 equal spheres and between 




FIG. 1: Geometry of the van der Waals interaction between 
2 spheres. 
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FIG. 2: Resonant plasmon frequencies as the function of nor- 
malized distance between two identical spheres (m = 0) with 
Drude's dispersion law. 



a sphere and perfectly conducting plane. 

Let us consider plasmon oscillations in a 2-sphere clus- 
ter with the geometry shown in Fig[TJ 

It is natural to use bispherical co-ordinates to calculate 
full plasmon spectrum in a two sphere cluster. The result 
of calculations is shown in Fig. [5] for axially symmetric 
modes (to = 0). Other modes have similar behaviour. 

From this figure one can see that there exist three 
different types of the modes: antisymmetrical (with re- 
spect to z) L modes and 2 symmetrical M and T modes. 
The L and T modes represent unbound states of plas- 
monic atoms in a single sphere while M modes represent 
the bound states of plasmonic atoms, that is plasmonic 
molecules These results only partially agree with 

previous calculations of plasmon spectrum [17|-[l9[ which 
do not predict plasmonic molecules at all. This difference 
is not very important in plasmon spectroscopy when a 2 
sphere cluster is excited by uniform (far) optical field 
with lj < uj p i j \/2, because plasmonic molecules do not 
exist under these conditions. However, when consider- 
ing van der Waals interactions ALL the modes should 
be taken into account , because the vacuum fluctuations 
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FIG. 3: Resonant plasmon frequencies as the function of nor- 
malized distance between two identical spheres (m=0) with 
Drude's dispersion law. 



are presented in ALL the modes. So one should expect 
substantial difference in van der Waals forces calculated 
within different approaches and (l7t-(l9|) 

Now consider the properties of different plasmonic 
modes in more details. The plasmon frequencies of L 
modes increase with the sphere separation and give rise 
to a van der Waals attraction. On the contrary, the plas- 
mon frequencies of M and T modes decrease when the 
distance between the spheres increases leading to a van 
der Waals repulsion. Symmetrical T modes are almost 
independent of the sphere separation and give no sub- 
stantial contribution to the van der Waals energy. Thus, 
L and M modes give the main contribution to the van 
der Waals forces. 

These features of plasmonic modes have a simple phys- 
ical interpretation (Fig. [3]). For example, for m=0 (axi- 
ally symmetrical modes), the dipole momenta of spheres 
for any mode is oriented along the z-axis. However, due 
to symmetry reasons, the dipole momenta of spheres with 
excited L modes are in the same direction, as shown in 
Fig. [3] , while the dipole momenta of spheres with ex- 
cited M or T modes have opposite orientation. Such ori- 
entations of dipole momenta result either in attraction 
between them, as is the case for the L modes, or in re- 
pulsion between them, such as for M and T modes. For 
the case m = 1 (angular dependence cosip or simp) we 
have similar sitiation. 

It is possible to find the following asymptotic expres- 
sions for the eigen- values of the dielectric constant in the 
case to ^ 1; L = 0, 1, 2, 3, . . . , mr) Q ~ 1: 



= - 



(L + m + 1) coth {{L + to + 1/2) %) + tanh r/ 
(L + to + 1) — tanh rj 

(2) 

where coshr; = Ri 2 / (2R ). 

An analogous situation takes place for symmetric M- 
modes. Asymptotic expressions for this case have the 
following form (to ^> 1,M = 0, 1, 2, ... , mrj ~ 1): 
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For a nanosphere and a perfectly conducting plate, 
plasmonic molecules can not be excited and only un- 
bound plasmonic atoms (anti-symmetrical L modes) will 
contribute to the van dcr Waals interaction 
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To find the van der Waals force between a sphere and a 
plate we derive ^ with respect to the distance between 
the sphere and its mirror image, that is R\2 in Fig. [1] 
where the x-axis should be thought to be replaced by 
the perfectly conducting plane. To estimate the van der 
Waals energy for small distances between the spheres, 
only large values of m are important and we can use the 
asymptotic solution ([2]) . Substituting ([2]) into ([4]) and 
subtracting the vacuum energy of free space ( R\2 — > oo), 
we obtain the following asymptotic expression 
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where A is the gap width between the sphere and its im- 
age. However, in the region where the sphere is not very 
close to the plane one has to calculate plasmonic frequen- 
cies of all L modes numerically and then to sum up the 
corresponding plasmon frequencies. The result of such 
summation for m = 0...50, L = 0...100 ( 2*5000 modes) 
is shown in FigfJ]( bottom lines). For the eigenvalue cal- 
culation we kept 1500 equations for each m = 0..20 and 
use asymptotics for m > 20. The dashed line corresponds 
to the asymptotic solution ([5]), while the solid line is the 
result of numerical calculation. The agreement between 
asymptote and numerical solution is very good for small 
distances. Our results are in agreement with the proxim- 
ity theorem but do not agree with results of recent paper 
[ 1 31 ] . where a more substantial singularity was predicted. 

The van der Waals interaction between two equal 
nanospheres is more complicated because now we have 
to take into account the zero-point energies of all modes, 
including plasmonic molecules (M modes) and unbound 
plasmonic atoms (L and T modes): 
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The contribution from L modes was already calculated 
in Eq. ([5]) , and we only have to find the contribution from 
M and T modes. Beforehand it is obvious that T modes 
will give only a small repulsive contribution to van der 
Waals forces C/ v T dw « Vf dW . For 
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FIG. 4: The van der Waals energy between 2 equal spheres 
and different contributions to it; Dashed lines correspond to 
asymptotic solutions ([5]) , (JSJ and their sum. 



only plasmon frequencies with high values of m will con- 
tribute to asymptotic expressions. We can, thus, use the 
asymptotic solution (|3|) to estimate the contribution from 
M modes to the van der Waals energy. Substituting © 
into ([7]), we obtain the asymptote 



UZ W « 0.08387^^ 



(8) 



The results from numerical summing up of the zero en- 
ergies of all M and T modes with m = . . . 50, M (T) = 
. . . 100 ( 4*5000 modes) and the corresponding asymp- 
totic expressions are shown in Fig [4] which shows a very 
interesting interplay between repulsive M modes and at- 
tractive L modes. Both sets of modes have the same 
singularity for close separation between spheres, but of 
opposite signs. The amplitudes of these singularities turn 
out to be almost equal in value (0.096 and 0.084). As a 
result, the full van der Waals energy between spheres 
remains attractive but small in comparison with the con- 
tributions from both the repulsive plasmonic molecules 
(M modes) and attractive unbound plasmonic atoms (L 
modes). It is interesting that the small contribution of 
symmetric plasmonic atoms (T modes) can not be ne- 
glected in this case. 

Our approach can be easily generalized to the case of 
unequal spherical nanoparticle of the radii R\ , i?2 . A 
complete analysis of this general case will be published 
elsewhere. Here we restrict ourselves to the case of closely 
spaced nanospheres. Again, the main contribution to 
the van der Waals forces is due to the modes with large 
m, and the contributions from plasmonic molecules and 
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FIG. 5: Dependence of functions /m and fi ([9| on ratio 
<^pi,i/^ P i,2 of the bulk plasmon frequencies 

unbound plasmonic atoms can be estimated as 

^vdw ~ nUpllJM [ Upl2 ) A (R 1+ R 2 ) (Q) 

U vdw ~ to pi ,2/i ^ — J A (i? 1+ fi 2 ) 

where and are the bulk plasmon frequencies 

of the spheres and the functions /m and Jl are shown in 
FigEl For equal bulk plasmon frequencies © reproduces 
(|5|) and ([8]). However, for different plasmon frequencies, 
the contribution from plasmonic molecules, which give a 
repulsive force, becomes less important in comparison to 
unbound plasmonic atoms, leading to an attractive force. 
This fact is in agreement with the case of a nanosphere 
near a perfectly conducting plane, where the plasmonic 
molecules do not contribute to the van der Waals energy 
at all. 

In the experiment our results could provide additional 
proof of the existence of plasmonic molecules by di- 
rect measuring of the van der Waals forces between the 
nanospheres, clarifying in this way an unsettled contro- 
versy with respect to their existence. Indeed, the cal- 
culations made in (l~7j-p^| do not predict symmetrical 



M modes in plasmon spectra. As a result, the full van 
der Waals energy will be described by the lower curve 
in FigdJ For example, for Al nanoparticles of 50 nm ra- 
dius with a gap of about A=lnm, van der Waals force 
within the plasmon spectrum [H-0] is about 1 nN while 
it is about 10 nN for the plasmon spectrum predicted in 
Such a difference is well measurable with AFM 
techniques. 

On the other hand, a substantial dependence of con- 
tributions to the van der Waals energy from plasmonic 
molecules and unbound plasmonic atoms on geometry 
and materials can have important impact in Nanotechn- 
logies. The relative contribution of repelling plasmonic 
molecules and attracting unbound plasmon atoms de- 
pends critically on the bulk plasmon frequencies of the 
nanoparticles. Therefore it would be possible to reduce 
substantially the van der Waals forces between different 
elements and at the same time, the unwanted phenomena 
(such as adherence and sticking of nanoparticles), if one 
provides the same plasmon frequencies for different ele- 
ments to enhance the contribution of repulsing plasmonic 
molecules. 

In conclusion we have calculated the energy of the van 
der Waals interaction between 2 closely placed metal- 
lic nanospheres as the energy of vacuum fluctations of 
ALL plasmonic modes existing in such system. The re- 
sults obtained depend crucially on the existence of re- 
cently predicted bound states of plasmons ("plasmonic 
molecules"). This allows us to suggest experimental mea- 
suring of van der Waals forces between 2 nanospheres as 
the proof of existence of plasmonic molecules. 

Acknowledgments 

The authors thank the Russian Basic Research Foun- 
dation (V.K., grants #07-02-01328, #05-02-19647), Pre- 
sidium of Russian Academy of Sciences ( V.K., Program 
"Quantum Macrophysics" ) , and the University Paris 6 
for partial financial support of this work. A.L. acknowl- 
edges partial financial support by the European Contract 
STRP-12142 NANOCASE. 



[1] V.V. Klimov and D.V. Guzatov, Phys. Rev. B 75, 024303 
(2007) 

[2] V.V. Klimov and D.V. Guzatov, Quantum electronics 37, 
209 (2007) 

[3] Y.C. Cao, R. Jin, and C.A. Mirkin, Science 297, 1536 
(2002). 

[4] D.J. Maxwell, J.R. Taylor, and S.M. Nie, J. Am. Chem. 

Soc. 124, 9696 (2002). 
[5] L. Martin-Moreno, F.J. Garcia- Vidal, H.J. Lezec, K.M. 

Pellerin, T. Thio, J.B. Pendry, and T.W. Ebbesen, Phys. 

Rev. Lett. 86, 1114 (2001). 
[6] E. Altewischer, M. van Exter, and J. Woerdman, Nature 

418, 304 (2002). 
[7] A.W.C. Lau, D. Levine, and P. Pincus, Phys. Rev. Lett. 

84, 4116 (2000). 
[8] C. Noguez, J. Phys. Chem. C 111, 3806 (2007). 
[9] H. Xu and M. Kali, Phys. Rev. Lett. 89, 246802 (2002). 



[10] P. Johansson and P.Apell, Phys. Rev. B 56, 4159 (1997). 
[11] L. H. Ford, Phys. Rev. A, 58, 4279 (1998). 
[12] D. J. Mitchell and B. W. Ninham , J. Chem. Phys. 56, 
1117 (1972). 

[13] C. Noguez, C. E. Roman- Velazquez, R. Esquivel-Sirvent 
and C. Villarreal, Europhys. Lett., 67 (2), 191 (2004). 

[14] G. Barton, Rep. Prog. Phys. 42, 963 (1979). 

[15] F. Intravaia, A.Lambrecht, Phys. Rev. Lett. 94, 110404 
(2005). 

[16] C. Henkel, K. Joulain, J.P. Mulet, and J.J. Greffet, Phys. 
Rev. A 69, 023808 (2004). 

[17] R. Ruppin: J. Phys. Soc. Japan 58, 1446 (1989) 

[18] P. Nordlander, C. Oubre, E. Prodan, K. Li, M.I. Stock- 
man: Nano Lett. 4, 899 (2004) 

[19] H. Wang, D.W. Brandl, P. Nordlander, N.J. Halas: Acc. 
Chem. Res., DOI: 10.1021/ar0401045 (2006) 



